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The Coulomb interaction between massless Dirac fermions may induce dynamical chiral symmetry 
breaking by forming excitonic pairs in clean graphene, leading to semimetal-insulator transition. If 
the Dirac fermions have zero bare mass, an exact continuous chiral symmetry is dynamically broken 
and thus there are massless Goldstone excitons. If the Dirac fermions have a small bare mass, an 
■ approximate continuous chiral symmetry is dynamically broken and the resultant Goldstone type 

excitons become massive, which is analogous to what happens in QCD. In this paper, after solving 
Dyson-Schwinger gap equation in the presence of a small bare fermion mass, we found a remarkable 
reduction of the critical Coulomb interaction strength for excitonic pair formation and a strong 
enhancement of dynamical fermion mass. We then calculate the masses of Goldstone type excitons 
using the SVZ sum rule method and operator product expansion technique developed in QCD and 
find that the exciton masses are much larger than bare fermion mass but smaller than the width 
of dynamical fermion mass gap. We also study the spin susceptibilities and estimate the masses of 
non-Goldstone type excitons using the same tools. 
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I. INTRODUCTION 



Graphene is a recently fabricated two-dimensional electron system. It has attracted a great deal of research activities 
since it exhibits many interesting properties and also has remarkable potential technical applications [l], 0. Its 
fundamental low-energy degrees of freedom are two-dimensional massless Dirac fermions, which obey relativistic 
'— Dirac equation. Due to the special linear dispersion of Dirac fermions, the density of states vanishes linearly near the 
. ' touching points of conduction and valence bands. Therefore, graphene is classified as a semi-metal and the Coulomb 
interaction between massless Dirac fermions is unscreened. 
\q ' The semi-metal ground state of graphene is stable when the Coulomb interaction is weak. Its microscopic action 
Q\ \ respects an exact continuous chiral symmetry. Due to the masslessness of Dirac fermions, graphene exhibits highly 
unusual behaviors, such as minimum conductivity il, 2j], quantum Hall effect Q, and marginal Fermi liquid behavior 
[|J. However, when the Coulomb interaction is sufficiently strong, a finite mass gap will be dynamically generated 
Q\ • due to excitonic pairing instability [5l4l2|. As a consequence, the continuous chiral symmetry is broken and the 
graphene undergoes a quantum phase transition from semi-metal to excitonic insulator. The mechanism underlying 
this phase transition can be considered as a concrete realization of the non-perturbative phenomenon of dynamical 
chiral symmetry breaking (DCSB) that was originally proposed by Nambu and Jona-Lasinio in the context of particle 
\ physics [l3j . 

Once a finite fermion mass is generated dynamically, both the ground state and the low-energy elementary ex- 
r> | citations of graphene are fundamentally changed. The ground state becomes insulating, and the only low-energy 
excitations are massless Goldstone bosons induced by dynamical breaking of exact continuous chiral symmetry. These 
bosons are composed of quasi-particles (fermion) and quasi- holes (anti-fermion), called excitons, and dominate the 
low-energy behaviors of graphene. 

Note that the excitonic pairing triggered by Coulomb interaction is not the only way to generate a fermion mass gap. 
For example, the Kekule type distortion can open a small fermion gap [Du - tlql. M oreover, the spin-orbit interaction 
may generate a small fermion gap that has opposite sign at K and K' points [17] . which was proposed to be able to 
turn graphene into a topological insulator [17| . Other gap generating mechanisms (for instance, geometry confinement) 
are also possible (l8| . From an experimental point of view, a finite fermion mass gap is expected to be present if 
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there is an exponential suppression of certain observable quantities, such as specific heat and susceptibility, at low 
temperature. However, the graphene may display different behaviors if the fermion gap is generated by different 
mechanisms. To see this clearly, it is helpful to make a symmetry analysis since different gap generating mechanisms 
are usually associated with different symmetry breaking patterns. 

When there is no excitonic pairing instability, the fermion mass may be generated by other mechanisms such as 
Kekule distortion. In this case, the continuous chiral symmetry is broken explicitly, rather than dynamically, and there 
are no Goldstone type excitons. If the fermion mass is completely generated by excitonic pairing instability, then there 
are massless Goldstone type excitons. The third possibility is that more than one mechanisms are important. Before 
the Coulomb interaction is turned on, the Dirac fermions may already have a finite mass due to certain mechanism. 
If this bare mass is small, the system possesses an approximate continuous chiral symmetry. When the Coulomb 
interaction is turned on, the Dirac fermions can acquire further mass gap due to excitonic pair formation. Once this 
happens, the approximate continuous chiral symmetry is dynamically broken (l9| . As a consequence, the Goldstone 
type excitons have finite masses [l9| . 

We believe that the third possibility is of particular interests for two reasons. First, this possibility can happen in 
realistic graphene materials. For instance, when Coulomb interaction and Kekule distortion (or other mechanisms that 
can induce bare fermion mass) are both present, there may be dynamical breaking of approximate chiral symmetry. 
Second, the physical picture of this possibility is very similar to what happens in QCD. In QCD, it is well-known that 
the u and d quarks have small bare masses. When these quarks acquire dynamical mass due to the formation of chiral 
condensate, the approximate chiral symmetry is dynamically broken and there appear massive Goldstone bosons, 
which are identified as the tt mesons [19| . The massive Goldstone type excitons generated in graphene correspond to 
the massive 7r mesons in QCD. 

In this paper, we study the dynamical breaking of approximate chiral symmetry in graphene. We first assume a bare 
constant mass for the Dirac fermions, and then calculate the dynamical fermion mass by solving the corresponding 
Dyson-Schwinger (DS) gap equation. From the solutions, the chiral condensate can be easily obtained. We will show 
that a small bare fermion mass can greatly catalyze the formation of excitonic pairs. The critical Coulomb interaction 
strength for excitonic pairing instability is reduced to zero, while the dynamical fermion gap is significantly enhanced. 

The dynamical fermion gap breaks the approximate chiral symmetry, so there appear massive Goldstone type 
excitons. Now a natural problem is to estimate the masses of these excitons, which can help to understand the low- 
energy excitations in the dynamical symmetry breaking phase. The close analogy between the approximate DCSB 
behaviors in graphene and in QCD allows us to calculate the exciton mass using the tools developed in QCD. In 
particular, we will use the Shifman-Vainshtein-Zakharov (SVZ) sum rule method and the operator product expansion 
(OPE) technique. OPE was proposed independently by Wilson [2(| and Kadanoff [2l[. Wilson's motivation came 
from the urgency to study the problem of strong interactions associated with the nuclear force, while Kadanoff applied 
the OPE to understand critical phenomena in condensed matter physics. Nowadays, OPE is known to be a powerful 
tool of quantum field theory describing elementary particle physics [22j and condensed matter physics [23| . SVZ sum 
rule technique was developed with the aim to understand strong interaction 0, HH . Among its wide applications 
in hadronic physics, SVZ sum rule is particularly powerful in calculating the masses of mesons and baryons that are 
bound states of quarks and/or anti-quarks. These methods were also applied to condensed matter systems, including 
degenerate electron gas [26[ and cold atom system [27|. In this paper, we calculate the masses of excitons using these 
analytical tools. 

These Goldstone type excitons are associated with the onset of charge density wave (CDW) or spin density wave 
(SDW) orders. In addition to these excitations, there are also collective spin triplet excitons, which are not related 
to any symmetry breaking. The masses of these non-Goldstone type excitons will also be calculated by means of the 
same SVZ sum rule and OPE methods. We hope these exciton modes could be found in future experiments (such as 
inelastic neutron scattering) in clean graphene sheet. 

The rest of this paper is organized as follows. In SecJlTl we give the continuum model of two-dimensional Dirac 
fermions interacting through Coulomb potential. We study the DS gap equation in the presence of small bare fermion 
mass and calculate the chiral condensate using the dynamical fermion mass. In Sec lIIH we study the Goldstone type 
excitons and calculate their masses using SVZ sum rule and OPE techniques. In Sec lIVi the masses of non-Goldstone 
type spin excitons are calculated. We give a brief summary and conclusion in SecfVl 



II. DYSON-SCHWINGER EQUATION AND DYNAMICAL MASS GENERATION 

Since the seminal work of Nambu and Jona-Lasinio, DCSB has been investigated in the context of particle physics 
for nearly fifty years [l3l |28 33] . In particular, it is one of the most prominent features of QCD. Unfortunately, the 
structure of QCD is too complicated, so the problem of DCSB in QCD has not yet been solved satisfactorily, although 
there have been remarkable progress in its supersymmetric version [34]. 
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In order to gain insights into QCD, some theorists turn to QED3. Despite its simple structure, QED3 shares a 
number of salient features with QCD: asymptotic freedom (35[, DCSB [33, and confinement [H, Specifically, 
Appelquist et al. found that DCSB can take place when the fermion flavor is less than certain critical value, Nf < Nj in 
QED3 j32|. Besides its relevance to particle physics, QED3 of massless Dirac fermions also has important applications in 
condensed matter physics. Indeed, it is the low-energy effective field theory for a wide class of planar strongly correlated 
electron systems, especially high temperature cuprate superconductor [38|]. The DCSB in QED3 is interpreted as the 
formation of two-dimensional Heisenberg antiferromagnetic order, which is the ground state of undoped cuprates [38| . 

In the case of graphene, the Coulomb interaction between Dirac fermions plays an essential role since it is poorly 
screened. There is a close similarity between the low-energy continuum theory of graphene and QED3 . It is therefore 
not surprising that Coulomb interaction can lead to DCSB. However, there is also important difference: the Coulomb 
potential is non-relativistic and contains only the temporal component of gauge field. Furthermore, although the 
massless Dirac fermions are two-dimensional, the electromagnetic field propagates in three spatial dimensions. 

Following Ref.JgJ, we describe the Dirac fermion in graphene by reducing QED4 to (2+l)-dimension. The effective 
action has the form 



S = I dtd rip s (t,r)(ij d t - ivF~f l ®i — m s )ip s (t,r) 

~ / d^d^'dWrV s (^,r) 7 V s (^,r)^o(^-^^|r-r'|)Vv(^^I■')7V S '(^^r'), (i) 



where the Fermi velocity Vp — c/300 and subscripts s,s' = l,2(or f, J.) are spin indices. The spinor field is ip = 
(Ak, Bk, Bk> , Ak>) with sublattice indices A, B and momentum valley indices K, K' . The 7 matrices are defined as 

1 2 \ ( -en \ / -er 3 \ / 1 2 



1 1, )■'■ u j' 73 = U3 ;' 75 -vo-i 2 

which are used by Q and [HI, and obey the Clifford algebra {7^,71/} = 2g M „ with fi, v = 0,1,2. Both 73 and 75 
anti-commute with 7^. Note that different choices of the 7 matrices do not change the final results obtained in this 
section. The bare Coulomb potential Uo(t, r) between Dirac fermions is given by 

%0 = ^, (2) 
er 

where e is the dielectric constant. In graphene, it is convenient to define a fine structure constant as a = e 2 /ehvp. It 
takes different values when graphene is placed on different substrates. For graphene on substrate Si02, a w 0.8; for 
graphene in vacuum (suspended), a = 2.2 [l~2l [l8|. 

When the Dirac fermions are massless, m s — 0, the Lagrangian (1) has an exact continuous chiral symmetry, 
%jj — > e^ 73,5 ')/'. However, this symmetry is not respected by the mass term ijji/j. Besides, when a finite fermion mass 
is generated by formation of excitonic pairs due to Coulomb interaction, this symmetry will be dynamically broken. 
Using the DS equation approach, it is argued that a sufficiently strong long-range Coulomb interaction can induce a 
dynamical fermion mass, {iptf)) ^ 0, thus leading to quantum phase transition from semimetal to excitonic insulator 

iiae. 

When the Dirac fermions have small bare mass m s with s = 1,2, the Lagrangian respects only an approximate 
chiral symmetry. We will solve the corresponding gap equation in the presence of a small bare fermion mass and 
calculate the vacuum chiral condensate (vac^s^/v |vac), which is essential for the later computation of exciton mass in 
the next section. Due to the above definition of gamma matrices, the mass term i/)ip corresponds to Kekule distortion 
d, Ojj]. More generally, Kekule distortion is formulated by (ReS)ip^ + (ImS)'ipi^5tp Q, but it can be transformed into 
rwipip with m — yj (Re<5) 2 + (Imi5) 2 = |5| by absorbing the phase into the 4-spinor tj) [HI ]. We should also note that 
the physical origin of bare mass do not affect the solution of gap equation and the magnitude of the corresponding 
mass generation only depends on the magnitude of the bare mass. For simplicity, we are mainly interested in the case 
where mi and 7712 are nearly the same, i.e., \mx — 7712I <S m with m = (mi + m,2)/2. 



The DS equation for Dirac fermion propagator of ipi reads 11 1 

(270 (L) 2 V{P ° ~ fc °' P ~ k ^( fc °' k )^' ^ 

where the bare propagator is 



(4) 
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and the full fermion propagator has the form 

S(po,p) = s t , r. (5) 

Hereafter we use the unit Vp = 1. Following previous works [5, 6], we adopt the instantaneous approximation and 
replace the gap function Ai(p , p) by Ai(0, p). Then we have 



\ / \ a f d 2 k 

A 1 (p) = m 1 + - I ^^(P,k) 



AlW (6) 



v /k 2 + A 1 (k)2 



with the kernel function 



where 



J(P ' k) = lp-^' (7) 



/%) = 3—3 T37TT- (») 



2m i _j_ q 2 -4m{ 

q 



arctan ( 4^- 

V 2mi 



Nf is the number of spin flavor and will be taken to be 2 in this paper. We approximate the kernel by its asymptotic 
value at p <C k and p > fc, so that 

J(p,fc) = 0(p-fc)ffl+0(fc-p)ffl (9) 
Therefore, the gap equation Eq.([6]) is written in the form 



Ai(p) = «!! + - / dkJ(p,k) 
k Jo 



Vfe 2 + Ai(fc) 2 ' 



(10) 



where A is upper momentum cut-off with the order of the inverse of the lattice constant. Here we choose its value as 
lOeV. 

In the present work, we solve the integral equation numerically. From A(p), the chiral condensate (vacl'f/'sVv |vac) 
can be evaluated by its definition! 11] 



(vac | -0s Vv | vac) = — tr Hm(vac|T-!/v(a;)'i/>s(0)|vac) = — tr 

x— >0 
A 



dp d 2 p 



A s (p)pdp 

7T Jo y / p 2 +A s (p) 2 



2ir (2tt) 2 i - A s 

(11) 



The numerical results for chiral condensate (v&c\?p s ?p s i |vac) = 8 ss i{il}tp) and dynamical fermion mass at zero mo- 
mentum A(0) are present in Table 1 for the case of zero bare mass m — and in Table 2 for the case of small bare 
mass m = 1.00 x 10~ 7 A. Comparing these two tables, it is easy to see that the small bare mass of Dirac fermion leads 
to a substantial enhancement of chiral condensate. When interaction strength a = 5, the dynamical mass A(0) in 
Table 2 is nearly twenty times larger than the value of dynamical mass in Table 1, while the chiral condensate is about 
one hundred times larger. When there is no bare fermion mass, there is no excitonic pairing and the Dirac fermions 
remain massless for a < 2.4. The critical strength a c of the Coulomb interaction lies in the interval of (2.4, 2.5), which 
is a little larger than the result of [TTJ] . However, when the Dirac fermions have a small bare mass, excitonic pairing 
and DCSB can take place for any finite value of a. For graphene placed on SiC>2 with a « 0.8 [TH, EH, the dynamical 
fermion mass A(0) = 2.15 x 10~ 6 A comes mainly from the formation of excitonic pairs since it is much larger than 
the bare mass 1.00 x 10~ 7 A. 

In summary, the small bare fermion mass has two effects: it reduces the critical Coulomb interaction strength a c 
and enhances the magnitude of dynamical fermion mass A(0). These effects are important from both experimental 
and technical points of view. On the one hand, the fermion gap can be detected unambiguously in experiments only 
when it is sufficiently large [l8| . On the other hand, the graphene material with a large gap will have more technical 
advantages than that with a negligible gap fl8j ]. 

We would like to emphasize that, although the bare fermion mass can catalyse the generation of dynamical fermion 
mass, the latter is physically different from the bare mass. In this work, the bare mass is assumed to be generated by 
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several possible mechanisms and small in quantity. It can be non-zero even when the Coulomb interaction between 
Dirac fermions is completely ignored. However, the dynamical fermion mass originates from the formation of chiral 
vacuum condensate driven by Coulomb interaction and thus is a typical effect of strong correlation between fermions. 

It is also interesting to compare the excitonic pair formation with the Cooper pair formation in the BCS theory of 
superconductivity. Historically, the proposal of dynamical chiral symmetry breaking of Nambu and Jona-Lasinio was 
motivated by the BCS theory. However, there are essential differences between them. An excitonic pair is composed 
of a particle and a hole, and thus is neutral. On the contrary, a Cooper pair is composed of two electrons, and thus 
carries negative charge — 2e. Moreover, the formation of excitonic pairs breaks chiral symmetry and leads to insulating 
behavior, whereas the formation of Cooper pairs break local gauge symmetry and lead to superconductivity. 

The DS equation with a Kekule-distortion induced fermion mass was studied previously in Ref. [l6[ , which concludes 
that the dynamical mass generation due to interaction is independent of the homogeneous Kekule distortion. In 
Ref. [l6j . an important claim was that the gap equation is dominated by the large momentum regime. However, it is 
known that dynamical chiral symmetry breaking is a non-perturbative, low-energy phenomenon. This phenomenon 
can happen only when the interaction between fermions is weak at high energy/momentum regime and strong at 
low energy /momentum regime (i.e., asymptotic freedom). Therefore, the processes with small energy/momentum 
transfer should play dominant role in the formation of excitonic pairing. According to our numerical computation, 
the dynamical mass generation is indeed significantly affected by the presence of homogeneous Kekule distortion whose 
catalytic effect can be readily seen from the comparison of Table I and Table II. 

TABLE I: Numerical results of chiral condensate without bare mass. Here, the chiral condensate is p = — (tfiip)/A 2 . 



a 


2.4 


2.5 


5 


10 


oo 


P 





1.34(-31) 


1.84(-9) 


2.22(-7) 


4.57(-6) 


A(0)/A 





2.67(-21) 


1.70(-6) 


4.37(-5) 


3.45(-4) 



TABLE II: Numerical results of chiral condensate with baremass m = 10 A. 



a 





0.8 


1 


2 


2.2 


2.4 


2.5 


5 


10 


oo 


P 


3.20(-8) 


5.40(-8) 


5.91(-8) 


8.57(-8) 


9.17(-8) 


9.78(-8) 


1.01(-7) 


2.13(-7) 


6.37(-7) 


5.12(-6) 


A(0)/A 


1.00(-7) 


2.15(-6) 


3.00(-6) 


8.37(-6) 


9.62(-6) 


1.09(-5) 


1.16(-5) 


3.12(-5) 


8.19(-5) 


3.67(-4) 



In the calculation presented below, we will need the normal-ordered chiral condensate (vac| : ^(O)-0(O) : |vac) = 
(vac|V'(O)'0(O)|vac) — (f2|-0(O)V>(O)|f2). Quite different from the non-perturbative symmetry-broken vacuum state |vac), 
is the perturbative vacuum state that is chiral symmetric [401] . As a function of a, the value of (vac| : ^(0)^(0) : 
vac) will be specified in Sec lIIII 

III. SVZ SUM RULE ANALYSIS OF GOLDSTONE TYPE EXCITONS 

Since the Dirac fermions have a small bare mass, the Lagrangian of graphene respects an approximate continuous 
chiral symmetry. As emphasized by Weinberg (ljjj . when an approximate continuous symmetry is broken, the Gold- 
stone bosons are no longer massless. Instead, these bosons are massive. In the context of graphene, the Goldstone 
type excitons induced by dynamical breaking of approximate chiral symmetry have finite masses. In this section, we 
calculate the masses of these excitons and compare them with the dynamical fermion mass. 

An exciton is a boson composed of a Dirac particle and a Dirac hole (i.e., an anti- fermion in the terminology of 
particle physics). It can be described by a composite operator of spinor field and its conjugate. For graphene, when 
chiral condensate occurs (^ip) 0, the chiral symmetry U(4) is broken down to U(2) x U(2) and the number of broken 
generators is 8. The 8 corresponding Goldstone bosons can be described by V'73®cr M '0 and 7/^75 ® er M V> [H,|Ilj], where 
<7o is unit matrix. The modes of ^175 €5 are bond-density-waves that mix the K and K' points, similar to the 
Kekule distortion mode [H|,[3!|. "073 <8> cr^ip is related to the CDW or staggered SDW excitations ||. 

In this section, we will calculate the masses of the lowest-energy excitons corresponding to (j>a = "073 <g> (jQip and 
4>i = "073 ® &iip using the SVZ sum rule method, which is analogous to the procedure of computing pion mass in 
particle physics (25|. This procedure can also be applied to the excitons associated with ^75 ® a^ip and the results 
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are the same, so we will not discuss the case of 1/^75 ® cr^tp. The basic idea of SVZ sum rule is to compute one 
particular physical quantity in two different ways and then extract important information of some parameter (such 
as the mass of a Goldstone boson) by equating the expressions obtained by different ways. In the present problem, 
the physical quantities to be computed are the correlation functions of composite fields cf>o and 4n. To apply the SVZ 
sum rule technique, OPE method will be used to separate the perturbative and non-perturbative contributions to 
the correlation function. The perturbative contributions are included in the so-called Wilson's coefficients and can 
be calculated perturbatively. The non-perturbative contributions are embodied as chiral condensates, which can be 
obtained from experimental data or calculated by some non-perturbative methods. In our case, the condensates are 
calculated by means of DS equation method. 

We first consider the field <j>Q. Its correlation function is defined as 

U(q) =i { d 3 :re l9X (vac|T0o(x)4(O)|vac), (12) 



where |vac) is the non-perturbative vacuum state in the chiral symmetry breaking phase. Substituting the 4 compo- 
nents of i/i, we have 



Af(K)A s (K) + At(K')A s (K') - B+ (K)B S (K) - B+(K')B S (K') = N A - N B , (13) 



which corresponds to the CDW excitation. In order to get information of the excitation, the correlate function will be 
calculated in two different ways: phenomenologically and theoretically. On the phenomenological side, the correlation 
function is related to physical observables, while on the theoretical side the same function is expressed in terms of 
fundamental parameters such as m and a which are treated as known numbers. Equating the results obtained in 
these two ways, we can obtain the expression for exciton masses in terms of fundamental parameters. 

To perform the phenomenological computation, we insert a complete set of physical states, 1 = J2 n \ n iPn)dT n (n,p n \, 
between the two operators in Eq. (112p . Here, p n is the momentum of the intermediate state and dr n — d 3 p n S(p^ — 
M 2 )9(p'^ l ) I (2tt) 2 denotes the integration over phase-space [42]. The index n can take both discrete and continuous 
values, which implies that both discrete bound states and continuous states are included. After integrating out space 
and momentum coordinates, we obtain 

ImII(g) =]>><5(<z 2 - M 2 n )F n , (14) 

71 

where F n = \(n, q|</>o|vac)| 2 . M n is the mass of the intermediate state \n) or the static energy in its rest reference 
system. For later convenience, we now introduce a function II(g 2 ), defined by 

U(q) = q 2 U(q 2 ), (15) 

and then have 

Imf^ 2 ) = TrU a 5(q 2 - m 2 J + ^j\f n 5{q 2 - M 2 ) (16) 

with F n — M 2 f n . Here we have isolated the contribution of the lowest- lying state (fiQ. The mass of (po, m V( , gives the 
position of the corresponding resonance. In the context of particle physics, f Vo is called decay constant of the bound 
state. In the context of graphene, it is related to the strength of exciton resonance. The second term of the righthand 
side of Eq. (fT6"|) includes both discrete and continuous states. Using the dispersion relation 

Il(q 2 ) = - / ds Ai^, 17 

7T J s- q- -ie 



we obtain 



n( g 2 )= 2 U ° +]^A^, (is) 

< - q 2 fM„ 2 - q 2 

where the first term is the contribution of the lowest-lying state and the second term is the contribution of higher 
excited states. 

In order to identify the mass m Vo of the lowest-lying exciton state corresponding to ^73 V'j we need to perform the 
theoretical analysis starting from the same correlation function. The correlation function can be calculated by means 
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of OPE method, which is able to account for the non-perturbative effects due to chiral-symmetry broken vacuum. We 
first write down the following expression 

i J d 3 xe^T0 o (x)^(O) = C (q) + C 2A (q) : ^(0) : +C 2 . 2 (q) : ^2(0) ■ +C 4 (q) : ^#^(0) : +... (19) 

The Wilson's coefficients C n 's contain the perturbative contributions, while the operators multiplying C n 's contain the 
non-perturbative contributions. From dimension analysis, we know that the mass dimension of the left hand side of 
Eq.(|19p is one, so each term in the right hand side should have the same dimension. Since the operators appearing in 
the right hand side are all local operators with increasing mass dimensions, the mass dimension of C n should decrease 
with n. This indicates that C„'s contain increasing powers of 1/q 2 as n grows. After taking vacuum expectation 
value, the contributions of higher terms are suppressed by (: iptp :)/q 2 for large q 2 . In the present work, we are mainly 
interested in the region —q 2 » (: (/'(/' :)| in Eq.(l9) and thus can keep only the first three terms. Generally speaking, 
C„ can be obtained by sandwiching the two sides of Eg. (1191) with a pair of states and comparing the results from 
both sides, as shown by Braateen [27j. Here for Co, we sandwich Eq. flT9"|) by perturbative vacuum state |f2), which 
eliminates the contributions of higher terms. Therefore Co will be calculated by evaluating the Feynman diagram 
of a fermion loop. For C2,i and C 2y2 , following the works of [IH, H^j, we sandwich Eq. (fT9]) by non-perturbative 
vacuum state |vac) and apply Wick's theorem to the left-hand side. Then a pair of fermionic operators are contracted. 
After calculating tree- level Feynman diagram and matching the two sides, we will get C 2 ,i and C 2>2 . Furthermore, 
perturbative calculations for Wilson's coefficients are performed in the framework of 1/Nf expansion and only leading 

term Cn 's are kept in the present work. 

Sandwiching Eq. (|19[) by the physical vacuum state, one obtains 

1%) = tf< 0) (g) + ($%){: :} + C<%){: :>. (20) 

Similar to the relation between II and II, we introduce C« 's satisfying 

C^(q)=q 2 C^(q 2 ), (21) 

and then have 

fiV) = C< V) + C!$(q 2 )(: fcfr :} + C^% 2 )(: ^ 2 :). (22) 
since can be obtained from evaluating the imaginary part of a fermionic loop diagram, we obtain 

ImC(°V) = -^=% 2 - 4m?) + ~^=0(q 2 - Am 2 ). (23) 

C^l (i=l,2) can be obtained by directly computing the tree diagrams, with the expression 



^% 2 ) = TrP^^^ (24) 
3(q z — ra^y 3q z 



in the region Q 2 = —q 2 ^> m 2 . 

Equating the right hand sides of Eq. fTS)) and Eq. (j2"2")l and substituting C„'s, we obtain 

7... ,^ /„ = i r ds im^(£j + 4m iHif : ) + 4w2 o ( ;^t 2 :) + o(-L). (25) 



Q 2 + Ko Z frQ 2 + M 2 nj s + Qi-te' 3(Q 2 ) 2 3(Q 2 ) 2 K Q< 

In order to extract the information of the lowest-lying state, it is helpful to introduce the Borcl transformation 

Bm*=~ lim (JL)» (26) 

Q 2 =nj\/ 2 ^oo n! dQ 2 

which can suppress the contribution from more massive states on the left hand side of Eq. (|2~5|) . Now we can obtain 
the sum rule 

U e- m *° /M2 + tfn*- M2JM2 =\l ddmCWWe-"" - |J + O(JL), (27) 



8 



where 



X 



-2mi(: V'iV'i :) - 2m 2 (: ^2^2 :}• 



(28) 



The contribution from higher states is suppressed exponentially and seems to be unimportant. However, these higher 
states include continuous ones, so they can not be simply neglected. We assume the continuous states begin at 
2A(p = 0) with A(0) being the dynamical fermion mass, and subtract their contribution from both sides of Eq. (j2"7} 
which is an analogue to the quark-hadron duality in particle physics [22|, |42|, [44], SB} • The resultant expression is 



(2A) 2 



dslmfl(,s)e 



-s/M 2 



(2A) 



dsImC^°\s)e 



-s/M 2 



/n , 1 N 
— F2+°( — )' 



3M 2 



(29) 



In this expression, X contains the normal-ordered condensates which were not specified in Sec. II. Here we defined 
them as 



/,7,1V 1 f 2Ali0) Mp)pdp 
(vac | : ipiipi ■ | vac) = ' 



7T Jo yV + Ai(p) 2 ' 



(30) 



and (: ^2^2 : ) can De defined similarly. From the numerical results of A(p) obtained Sec. II, the value of (: iptp :) can 
be computed by numerical integration. Then we obtain the sum rule formula 



(31) 



where ... is the contribution from the discrete higher states and is neglected altogether. Differentiating the above 
expression with respect to 77 = 1/M 2 and introducing the ratio between the derivative and Eg. (1311) . we eventually 
obtain the sum rule for the mass of exciton ipo : 



2 _ j££ v?e-"*«du - fm^ ijvh + (1 ^ 2) 

Vo ~ £ e ~ u2r>du + -)v + (1^2)' 



(32) 



Besides, the sum rule formula for f Va is 

r m 2 77 



1 

2^ 



2Ai 



2mi 



e- a "d U + 4T " l(: f 1 ^ : V (lo2) 



(33) 



We next consider the correlation function of <pi = - 0i73'02 ~ ^.{Al- — B^,B±. Analogously, we define 

U 1 {q)=i J d 3 xe l9X (vac|T^)<rf(0)|vac), (34) 

which is the counterpart of the transverse component of staggered spin susceptibility investigated in the context of 
high temperature cuprate superconductors [461 ] . Similar OPE is applied to Hi (q) , and the results for the corresponding 
coefficients are 



ImC#V) 



1 



6{q 2 - Am 2 ) 



4-yV 

with m = (mi + w 2 )/2. The resulting sum rule formula after Borel transformation is 



f<pi e 



^IJM 2 



1 

2^ 



2A 



e -" 2 /M 2 du + 



2 m 



(3m 2 - mi)(: ipiipi :) 



3M 



n- 



l(l«2) 



Differentiate it with respect to 77 = 1 /M 2 then one gets the sum rule formula for the mass of exciton ip\ : 

2 = lit u 2 e~^du - [f (3m 2 - mQQ fogi :) + (1 g 2)] 
/ 2 ^e-« 2 ^ U +[f (3m 2 -m 1 )(:^i^i :)i/+(l«2)] 

with 2m = mi + m 2 and 2A = Ai(0) + A 2 (0). 



(35) 



(36) 



(37) 
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The sum rule formulae for the masses of excitons fa and fa can be derived similarly. The formula for m ip2 is exactly 
the same as m lpi , Equ. (|37[ ). For fa, one considers the following two-point correlation function 

113(9) = * J dW 9 *(vac|T0 3 (x)0 3 (O)|vac), (38) 

where fa is the longitudinal component of the staggered SDW, given by 

fa~(N At -N M )-(N Bf -N Bi ). (39) 

The sum rule formula for m Va is the same as m Vo , which is presented above. 

In our truncated computation, these sum rule formulae for the exciton masses Eq. (|32[) and Eq. (|37p arc functions of 
the Borel parameter M or 77. In order to fix the value m v for exciton tp, we need a criterion for the choice of r\ and 
the corresponding m v (rf). In the sum rule literature, the optimum estimate- value may be determined by an extreme 
value, an inflection point [H, S3 or a plateau [42| of the function. In the present work, we choose a plateau near the 
extreme value of the function to be the final output for m v or f v . 

For mi = rri2 = 1.00 x 10~ 7 A and a = 0.8 (SiC>2 case), the numerical results for the chiral condensate and mass gap 
are -(: ■) = -{■ fa^2 :> = 1-70 x 10~ 12 A 2 and A(0) = 2.15 x 10~ 6 A. The cut-off scale A is normally set as lOeV. 

From these quantities, the SVZ estimates for the masses of the excitons can be made with the criterion discussed 
above: ipo, ip\, (f2 and ips have the same mass 0.028meV. For mi ^ 7712, the four- fold mass degenerate is lifted to 
two-fold: ip 1 and (p2 have the same mass, which is a little smaller than and m (p3 . For the case of suspended 
graphene, a = 2.2, we have — (: ■tpi/j :) = 3.33 x 10~ n A 2 , A(0) = 9.62 x 10~ 6 A. Now the mass of ip particles is about 
0.098meV. 

It is now necessary to summarize and compare the relevant energy scales discussed above. The fundamental energy 
scale in the present problem is the ultra-violet cut-off A = lOeV in graphene, which is normally determined by the 
lattice constant. In this paper, we assume a bare fermion mass m = 1.00 x 10~ 7 A, which may be generated by Kekule 
distortion or other mechanisms. All other energy scales are derived by explicit calculations. We note here that the 
exciton mass is smaller than 2A(0), which implies that the binding energy of the bound state is negative. 



IV. SPIN SUSCEPTIBILITY AND NON-GOLDSTONE TYPE EXCITONS 



In addition to the CDW and staggered SDW excitons studied in the last section, there is another kind of low-energy 
collective excitations: spin excitons. These spin excitons are not generated due to chiral symmetry breaking and thus 
are non-Goldstone type bosons. In the absence of any (bare or dynamical) fermion mass, these spin excitons are 
massless in clean graphene (4|| Sil and become massive in doped graphene • In the present problem, the fermions 
have finite dynamical mass in the chiral-symmetry broken phase, so the spin excitons are also massive. In order to 
study these triplet spin-1 excitons, we turn to study the spin susceptibilities which may be measured by inelastic 
neutron scatterin g [49j. These quantities were studied previously in an effective QED3 theory of high-temperature 
superconductors [51|. 

In this section, we will calculate the masses of spin excitons using the methods presented in the last section. To do 
this, we first define the spin operator as Si = ^t^ab^b, where i = x,y, z and a % are Pauli matrices. The transverse 
and longitudinal spin susceptibilities are defined as 

X+-(q)=% J d 3 xe iqx {vB,c\TS + (x)S-{0)\vac), (40) 

and 

Xzz(q) = i J d 3 xe iqx (v&c\TS z (x)S z (0) | vac) , (41) 

respectively. Here, |vac) is the non-perturbative physical vacuum state of the system in the chiral-symmetry broken 
phase. The quantity x h can be similarly studied and will not be discussed here. 

The spin susceptibility will be computed by means of SVZ sum rule method. We note that Rosenfelder applied the 
SVZ sum rule method to analyze density-density correlation function and estimate the position of plasmon resonance 
in a non-relativistic electron system [26| . In our case, the spin susceptibility will be analyzed in both phenomenological 
and theoretical ways, analogous to what we have done in the last section. 

On the phenomenological side, following the previous work (48l - l5l1 |. we assume that there is a resonance below the 
spin gap. This assumption is embodied in the expression 

-Im Xl (p) = FJ(p 2 - m 2 ) + V F n S(p 2 - M 2 ) + 0(p 2 - t/ sg W(p 2 ), (42) 

7T * ' 

n 
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where the first ^-function type resonance a corresponds to the so-called spin exciton. a = +, — , z for I = + — , — h, zz 
respectively. The spin gap U sg is the lower bound of the continuous spectrum. The second term is the contribution 
from the discrete excited bound-states, whose masses ^i n 's are either between m a and U sg or beyond U sg . From the 
dispersion relation formula, we obtain that, for p 2 < 0, 

Xi{P = - / ds - 

T Jo s-p z -le 



m 2 _ p 2 ^ M , _ p , j g _ p 



To compute the same correlation function theoretically, OPE is adopted to take into account the non-perturbative 
effects due to chiral symmetry breaking. For transverse spin susceptibility, we have 

i f d 3 xe i r x TS+(x)S-(0) = D (p) + D 2A (p) : &^(0) : +D 2>2 (p) : ^(0) : +D 4 (p) : ^^r>(0) : +... (44) 



As explained in Sec. Ill, the operators in the right hand side are all local operators, so the contributions of this series 
are suppressed by powers of (: ijjip '.)/q 2 . Therefore, we can keep only the first three terms. Sandwiching Eq. (|44l) by 
the physical vacuum state, we have 

X+ -(p) = D ( °\p) + D { 2 %)(: V^i :) + I>g(p)<: :> + ... (45) 

D n can be obtained by perturbative expansion computation and each D n is a power series in coupling a, therefore 
only the leading terms of are important. For Dq , its imaginary part is 

ImD'V.P 2 ) = e( P 2 - 4m 2 )-Pj= + Q(p 2 - V)-^. (46) 



8 v p 2 2p 2 v p 2 

After calculating the coefficients -D„'s in Eq. (|4"5]) . we get the truncated expression of transverse spin susceptibility for 
the case ra\ = 7712, as follows 



„ 2 _ 1 r dc ImD^(s) + 8m(: ^ :) p 2 
n J s — p 2 — ie 3(— p 2 ) 2 



where (: t/;^ :} = (: V'lV'i : ) = ( : 1 /'2'02 '■) and — p 2 3> tj 2 . Note that all terms with higher order of m 2 /p 2 are ignored. 
Equating the right hand sides of Eq. (|4"3|) and Eq. (|47|) , we get a sum rule formula 



1 [°° Im X+ _( g ) 1 f°° j IpxPf(s) , 8m(:#:)p 2 

ds 5 = - / ds + — 2T2 — • ( 48 

s — p 2 — ie it J s — p 2 — ie iy—p^y 



7T 



In order to subtract the contribution from the states beyond the spin gap in both sides of the above equation, duality 
is used and the upper cutoff is replaced by (2A) 2 , so that 

1 r 2 ^ Imx r (s) = - r ? ^ Im 1 0)(s) + M ;*t; 2 )pa . m 

7r J s — p 2 — ie 7r J s — p 2 — ie 3(— p 2 ) 2 

In this expression, the value of (: "ijiip :) was given by Eq. (|30[) . Now we can obtain the following sum rule formula 

F+ _ 1 f^ 2 Im^ 0) ( S ) 8mp 2 1 A(s)ds 



m+ -p 2 7T J s-p 2 - ie 3(-p 2 ) 2 it J ^ s + A(s) 

where ... stands for the contribution from the discrete excited bound-states. Furthermore, in order to extract the 
information of the lowest-lying state, Borel transformation can be used to suppress the contribution from excited 
states in the left hand side of Equ. ([5U)l . leading to 
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FIG. 1: Typical mass scales for effective fermion gap 2 A, CDW exciton m v and spin exciton m CT . The bare fermion mass m is 
chosen as 1/ieV. 



The second term of the left hand side of the above equation is the contribution from the more massive states, which is 
suppressed exponentially and will be neglected altogether in the following. Differentiating the above expression with 
respect to r\ = M 2 , we get the following sum rule formula for mass m + of spin exciton: 

The sum rule formula for m_ is exactly the same as m + . In addition, the longitudinal spin susceptibility \zz can 
be analyzed similarly. When — p 2 3> m 2 , the OPE result of Xzz(p 2 ,p) is 

X.,tf,V) = ~ r^ ImD f (S) + 8mi(: ^- )+ 9 8 r (: ^ 2:) P 2 - (53) 



tt Jo s — p 2 — ie 3(— p 2 

The final sum rule formula for the associated quantity m 2 will have the same form as Eq. (|52p in the case of mi = mi- 
Apparently, the three spin excitons are mass-degenerated. 

In our trancated computation, the sum rule formula for the exciton masses Eq. (1521) is function of the Borcl parameter 
M or rj. These two parameters are determined in the same way presented in Sec. III. For mi = 7712 = 10" 7 A and 
a = 0.8 (SiC>2 case), the numerical results for the chiral condensate are — (: ipiipi :) = — (: ^2^2 '■) — 1-70 x 10 -12 A 2 . 
From these quantities, the SVZ estimation for the mass of the spin excitons is m a = 0.034meV. For suspended 
graphene with a — 2.2, — (: ipip :} = 3.33 x 10~ n A 2 and the estimated value for the spin exciton mass is 0.12meV. 

The typical masses obtained above are illustrated in Fig.l. Comparing these masses, we see that the spin excitons 
are more massive than Goldstone type excitons, while their masses are smaller than 2A(0), which is consistent with 
the assumption made before Eq. p2")) . It is interesting to compare our results shown in Fig.l with the Fig. 3 of Ref. 
[i"o| , where the masses of pion-like excitons were calculated in the framework of lattice gauge theory. 
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V. SUMMARY AND CONCLUSION 

In this paper, we studied the dynamical breaking of approximate chiral symmetry in graphene and calculated the 
spectra of both Goldstone type and non-Goldstone type excitons. In the presence of a small bare fermion mass, the 
critical Coulomb interaction strength a c for dynamical fermion mass generation is reduced and the dynamical fermion 
mass is substantially enhanced. In particular, we found that the effective fermion mass (A(0) = 2.15 x 10 _6 A for 
graphene on Si02 with a = 0.8 and A(0) — 9.62 x 10~ 6 A for suspended graphene with a = 2.2) is much larger than 
bare fermion mass (to = 10 _7 A). Apparently, the enhancement of fermion mass originates from the excitonic pairing 
instability due to Coulomb interaction. 

When the fcrmions have small bare mass, an approximate chiral symmetry is dynamically broken, thus there appear 
massive Goldstone excitons. In the symmetry-breaking phase, there are no massless excitations, and all fermionic and 
bosonic excitations become massive. We calculated the masses of Goldstone type excitons using the SVZ sum rule 
method developed and widely used in QCD and show that they are larger than bare fermion mass but smaller than 
dynamical fermionic gap 2A(0). In order to take into account the effect of chiral symmetry breaking, OPE technique 
was used in the calculation of two-point correlation functions. In graphene, besides Goldstone type excitons, it is also 
interesting to study the non-Goldstone type, spin excitons. We specified their positions by the same SVZ sum rule 
method and found that the masses of these spin excitons are much larger than bare fermion mass but smaller than 
2A(0). Moreover, their masses are larger than those of Goldstone type excitons. 

In the theoretical treatment of graphene, the most widely used approximation is to keep only the nearest hopping 
and expand the fermion energy around the neutral Dirac points. Additional terms will be included in the effective 
continuous field theory when higher order corrections are taken into account. Some of these terms may break the 
chiral symmetry. Generally, the symmetry-breaking terms can appear in two classes: either as quadratic terms ifjTip 
or as quartic terms (ipTip) 2 with T = 1, 73, 175. The former corresponds to fermion mass terms and through a unitary 
transformation of the 4-spinor ip the coefficients of ij)Ttjj can be absorbed into the bare fermion mass m s introduced 
in our Eq.(l). Since our calculations and results depend only on the magnitude of to s , not on its physical origins, 
such terms will not qualitatively affect our conclusion if they are sufficiently small. The quartic term (4>ip) 2 can 
be considered as a Hubbard-type short-range interaction term. Its effect is more complex than the quadratic term 
because we need to study the interplay of the long-range Coulomb interaction and this short-range interaction pj. 
When such short-range interaction is weak, its contribution to the Dyson-Schwinger gap equation (3) can be studied 
by the methods presented in [7]. The calculations of the masses of Goldstone bosons and non-Goldstone bosons can 
be performed using the sum rule and OPE methods of Sec. Ill and Sec. IV, because these calculations rely only on 
the chiral condensate, which is easily obtained from the fermion gap function. However, if this quartic interaction is 
strong, the whole Lagrangian of the continuous theory is no longer chiral symmetric, and there will not be dynamical 
chiral symmetry breaking and Goldstone bosons. Indeed, our calculations are based on the assumption that the 
symmetry-breaking interaction is absent or sufficiently weak. 

After obtaining the masses of various types of exciton, the next problem is to judge whether and when these 
collective modes exist in graphene. Since their masses correspond to the resonance positions in the low energy region, 
we hope that NMR and neutron scattering might be able to address this problem, similar to the efforts in high 
temperature superconductors (5l| - [53l |. To make a connection with experiments, it is also necessary to calculate some 
observable quantities that can describe the effects of massive excitons. This issue is beyond the scope of the present 
paper and will be discussed in the future. 
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